"get.QSR" = function(sim.values, lag.vector){
		statistic.vector = NULL;
		pval.vector = NULL;		
		n=length(sim.values);
		#get autodistance-covariance at specific lag
		adcov=auto.dcov(sim.values,sim.values,max(lag.vector));
		V2 = adcov$adcov^2;
		S2 = adcov$S2;
		ns = adcov$ns;
		for (i in 1:length(lag.vector)){
			statistic.temp = ""
			for(k in 1:lag.vector[i]){	
				test.stat = ns[k]*V2[k]/S2[k]
				statistic.temp = paste(statistic.temp,",",test.stat,sep="");
			}
			statistic.vector = c(statistic.vector,statistic.temp);
			#other vectors  unavailable
			pval.vector = rep(-1,length(lag.vector));
			df.vector = rep(-1,length(lag.vector));
		}
	return(structure(list(name = "QSR",STATISTIC = statistic.vector,lags=lag.vector,PVAL=pval.vector,df=df.vector)));
}